m1a <- lm(log1p(abs(stan_violence$days)) ~ age+pid2+educ4+gender, stan_violence, weights = weight)
summary(m1a)

xtable::xtable(car::Anova(m1a))